Local modes, phonons, and mass transport in solid 4 He. 
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We propose a model to treat the local motion of atoms in solid 4 He as a local mode. In this model, 
the solid is assumed to be described by the Self Consistent Harmonic approximation, combined with 
an array of local modes. We show that in the bcc phase the atomic local motion is highly directional 
and correlated, while in the hep phase there is no such correlation. The correlated motion in the 
bcc phase leads to a strong hybridization of the local modes with the Ti(llO) phonon branch, 
which becomes much softer than that obtained through a Self Consistent Harmonic calculation, 
in agreement with experiment. In addition we predict a high energy excitation branch which is 
important for self-diffusion. Both the hybridization and the presence of a high energy branch are a 
consequence of the correlation, and appear only in the bcc phase. We suggest that the local modes 
can play the role in mass transport usually attributed to point defects (vacancies). Our approach 
offers a more overall consistent picture than obtained using vacancies as the predominant point 
defect. In particular, we show that our approach resolves the long standing controversy regarding 
the contribution of point defects to the specific heat of solid 4 He. 



PACS: 67.80-s,67.80.Cx,67.80.Mg 
I. INTRODUCTION 

Atomic motion in solid 4 He was extensively studied 
over the years, both theoretically and experimentallyfil. 
The quantum nature of the solid comes into play mainly 
through the large amplitude of zero point motion of the 
atoms in the lattice, an effect especially evident in the 
solid at its lowest density. In this paper, we shall discuss 
only these low density solids (~ 21 cm 3 molar volume). 
Calculations of the phonon spectrum within the self con- 
sistent harmonic approximation seem to agree fairly well 
with experiment in the hep phaseel, while in the bcc 
phase, only the energy of the Ti(llO) phonon is over- 
estimated by a factor of twctl Corrections resulting from 
the local motion of the atoms due to quantum-effects are 
treated using correlated basis wave functions!! A prob- 
lem which so far was not addresed within these models is 
that of mass transport, arising from atomic exchange at 
low temperature. At higher temperatures, this exchange 
is assumed to occur with thermally activated point de- 
fects, traditionally taken as vacancies. With point de- 
fects, some additional complexity is to be expected, as 
here one deals with objects of a size similar to that of 
the zero point amplitude of the atom. The existence 
of thermally activated non-phonon excitations is well es- 
tablished exnexjmentally by several-different techniques, 
such as NMRd&eIl , mass diffusionErO, ion mobilitvE9, and 
X-ray diffraction!! However, the properties of point de- 
fects deduced from these various experiments have so 
far led to controversial conclusions. For example, taking 
the experimentally determined activation energies at face 
value, the calculated density of vacancies is so large that 
their contribution to the specific heat should be compa- 
rable to that of the phonona!3fl This is not borne out 



by the specific heat data. Several attempts were made to 
reconcile this controversy by using the assumption that 
vacancies ia rSftUd 4 He are delocalized and occupy an en- 
ergy bancLj'Efra. In this approach, the magnitude of the 
specific heat is divided by the bandwidth, so that it can 
be made small through the choice of a suffficiently large 
bandwidth. However, in order to obtain a realistic magni- 
tude of the specific heat one has to chose an unphysically 
large bandwidth, so thai, this approach is unsuccessful in 
settling the controversyll 

In this paper, we propose a somewhat different de- 
scription of the solid, one which treats the motion of the 
atom inside it's potential- well as a "local mode". Lo- 
cal modes as excitations in anharmonic crystals were dis- 
cussed previously by Sievers and Takenotj, and by Flynn 
and StonehamEj. We show below that the other degrees 
of freedom are largely decoupled from this local motion. 
As such, this mode can be assumed to be an indepen- 
dent degree of freedom of the solid. The excited states of 
this local mode can be thermally activated, and seems to 
reproduce all the physical effects usually attributed to va- 
cancies. Moreover, we show that this picture has none of 
the inconsistencies connected with vacancies. The prop- 
erties of the local modes are different in the bcc and 
the hep solid phases. In the bcc phase, we show that 
it is energetically advantegous for the local motion of the 
different atoms to be correlated. The new ground-state 
has several new features, among them a softened trans- 
verse phonon mode and a high-energy excitation branch 
which opens an additional channel for atomic diffusion. 
In contrast, the triangular symmetry of the hep lattice 
frustrates the possibility of long-range ordering of the lo- 
cal modes, which therefore remain uncorrelated. In both 
phases the local modes are highly directional due to the 
anisotropic potential seen by an atom in the crystal, a 
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constraint which we show leads to a small contribution 
to the specific heat of the solid. 



II. LOCAL MODES IN BCC HE 

Calculations of the ground-state energy of bcc 4 He usu- 
ally employ variational wavefunctions that account for 
the short-range correlations between the atomscl. Since 
we want to focus here on the nature of the local motion 
of the atom inside the lattice, we will take a different 
approach. In this approach we would like to isolate the 
lowest energy excited state of the atom inside its poten- 
tial well, and treat it as a local excitation of the lattice. 
This excited state consists of a local oscillatory motion 
of the atom along a particular direction and produces 
an oscillating electric dipole, similarly to the usual Van- 
der Waals interaction. Unlike in the case of the Van-der 
Waals interaction, in which the dipolar fluctuations are 
random, we show that in the bcc solid these local dipoles 
are correlated and a new ground-state of lower energy 
is created. We therefore begin by investigating the po- 
tential well of an atom in the bcc lattice. This can be 
done in the followipg way: Using the standard helium 
pair-potential v(rp and taking the atoms as stationary 
we can map the potential-well near an atom along any 
direction in the lattice. In Fig.l we plot this potential 
along the main directions (100), (110) and (111). It is 
clear that in the directions normal to the cube's faces 
(i.e. (100), (010) etc.) the confining potential well is very 
wide with a pronounced double-minimum structure. 

We also plot in Fig. 2 the lowest two energy levels of a 
one-dimensional Schrodinger's equation for a 4 He atom 
solved in each of the potential wells. The energy differ- 
ence between these two levels is lowest in the (100) direc- 
tion, of the order of 10K. The atomic displacement is de- 
scribed by a mixing between the lowest two energy levels 
in the potential well, which corresponds to a motion with 
amplitude of ~ 1A (in the (100) direction). Although 
the above treatment assumes that other atoms are sta- 
tionary, it indicates the directions along which there will 
be large amplitude local motion in the solid. It is evi- 
dent that the motion will be highly directional, with the 
largest amplitude along the (100) direction, while in the 
other directions the amplitude is much smaller. Allow- 
ing other atoms to move as well, so that they get out of 
each other's way, can only soften the potential well and 
reduce the value of the energy level difference, since the 
potential will effectively become shallower. For exam- 
ple, a calculation which allows opposite atoms to move 
slightly, reduces the energy difference in the (100) direc- 
tion to almost 6K. 

In order to obtain meaningful numerical values, we 
scaled the potential v(r) so that the calculated one- 
dimensional potential energy and kinetic energy, aver- 
aged over the different directions, reproduce the known 
kinetic, potential and total energy of the bcc phasdil: 



(K.E.) ~ 34K, (V) ~ -40K, (E total ) ~ -6K. 

Based on the above calculation, we shall assume that 
the atoms have a local-mode that is highly directional 
along one of the directions equivalent to (100), and of 
energy 6 — 10K. Experimental evidence for the existence 
of such a "local mode" comes from NMR measurements 
of the linewidth (I/T2) in bcc 3 He- 4 He mixture crystalsO. 
The motional narrowing of the NMR resonance line with 
temperature shows thermall;v- 1| activated behaviour with 
an activation energy of TilKOu. At temperatures above 
IK, the NMR line in the solid becomes narrower than 
that of the liquid, indicating that atomic motion in the 
solid is faster than in the liquid. At the same time, the 
diffusion coefficient of the solid remains several orders 
of magnitude smaller than that of the liquida, indicating 
that this rapid motion is of a local nature. We propose to 
identify this rapid motion as associated with the excited 
state of an atom in the well, with an activation energy of 
7±1K. 

The highly directional motion of the atom leads to the 
creation of a local (oscillating) electric dipole in the di- 
rection of the motion. This local electric-dipole is created 
due to inertia, as the electronic cloud can be thought of as 
being slightly displaced relative to the ion. This electric 
dipole due to mixing of the \s) and \p) electronic-levels of 
the He atom. The amount of mixing can be estimated 
from perturbation theory as 

V; = \s) + \\p)^E ~ (V \E\ i&> - (a \E\ s) ~ A 2 (p \E\p) 
=> A 2 ~ 7/2.46 • 10 4 ~ 0.00284, A ~ 0.0168 (1) 

where \s) and \p) stand for the ground-state and first 
excited-state of the 4 He atom, A is the mixing coefficient 
and (p\E\p) ~ 2.46 ■ 10 4 Kjs the excitation energy of 
the first atomic excited-statetLS. The estimated mixing is 
small and the magnitude of the induced dipole moment 
is therefore 



H = e (V> \x\ ip) ~ 2eA (s \x\ p) ~ e • 0.03A 



(2) 



where (s\x\p) ~ 0.9A . The estimation of the mix- 
ing A and the dipole-moment |//| serves to set an upper 
bound on the magnitude of this effect. Since the atoms 
possess an oscillatory electric-dipole moment they have 
long-range dipole-dipole interactions. The instantaneous 
dipolar interaction energy is given by 
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dipole 



E 



3 cos 2 (p, ■ (r - r;)) - 1 



(3) 



where the sum is over all the atoms in the lattice, 
being the instantaneous coordinate of the z-th atom. For 
oscillating dipoles with random phases, the average in- 
stantaneous interaction energy summed over the lattice 
would be zero. However, the energy of a dipolar array 
can be made lower by correlating the phases of the oscil- 
lating atoms. The lowest interaction energy arrangement 
of the dipoles in the bcc lattice is such that they are os- 
cillating with the same phase. Since the direction of the 
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dipole shows the instantaneous direction of the motion 
or displacement, such a state is just a uniform motion 
or translation of the entire lattice. This arrangement is 
therefore unphysical, and we have to look for symmet- 
ric arrangements with respect to the number of up/down 
dipoles. The two arrangements shown in Fig. 3 are the 
two 'antiferroelectric' configurations along the symmetry 
axes of the crystal with individual dipoles oriented along 
the (001) direction, and a zero total dipole moment. For 
these two physical possibilities the sum in (^) with a unit 
dipole is given in Fig. 3. Thus, the ground state in our 
picture has the atoms executing this local oscillation in 
a correlated fashion, as shown in Fig. 3a along one par- 
ticular direction. Similar correlated motion exists in the 
other two orthogonal directions. 



III. ELEMENTARY EXCITATIONS OF THE 
DIPOLE GROUND-STATE 

The ground state of the dipoles described in the pre- 
ceding section will be affected by the excitations of the 
lattice, namely phonons. In fact, our basic assumption 
in which the local motion can be separated from these 
other degrees of freedom needs justification. The oscilla- 
tory atomic motion induced by the phonons will modu- 
late the relative phases of the dipoles. Let us look at the 
ground state of the dipoles concentrating on oscillations 
oriented along the (001) direction. We now need to con- 
sider only phonons which will modulate the local motion 
along this direction. In the bcc structure, only 3 phonons 
fulfill this condition: L(001), T(100) and Ti(110). Let us 
calculate the energy of the dipolar array when modulated 
along these 3 directions. For a modulation along so 



direction k , the dipolar interaction energy is given b; 



X(k) 



3 cos 2 (n ■ (r - r;)) - 1 



exp [2nik ■ (r 
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(4) 



At k = the interaction matrix X{k) is just the dipolar- 
energy (|). 

In Fig. 4 we plot the value of X(k), the energy of the 
dipolar array modulated by the relevant phonons: L(001), 
T(100), and Ti(110), for dipole moment = 1. We see 
that for a modulation by L(001) and T(100) the peri- 
odicity of X(k) is over a full unit-cell, that is twice the 
periodicity of these phonons. Since symmetric functions 
of periodicities Tr/a and 2tt / 'a are orthogonal, the wave- 
functions of the phonons and of the dipole-excitations 
are orthogonal along these directions. The dipole array 
cannot therefore be excited by any phonon along these 
two directions. Regarding the modulation by the Ti(110) 
mode, here the periodicity of X(k) is the same as that 
of the Ti(110) phonon, which can therefore couple to 
the dipole array. We conlude therefore that the coupling 
of the local modes to the lattice excitations is limited 



to a single phonon mode, justifying our assumption that 
the local modes can be treated separately to a good ap- 
proximation. Thus, the only elementary excitations of 
the dipole array would be in the (110) direction, in the 
form of the Ti(110) phonon. We shall now calculate the 
dispersion relation of such an excitation by a mean-field 
solution of an effective Hamiltonian. 

The Hamiltonian treatment of interacting local excita- 
tions was developed originally by HopfieldEJ for the prob- 
lem of excitons in a dielectric material. The local excita- 
tions are treated as bosons aad the effective Hamiltonian 
describing their behavior ial3 



(5) 



Hioc = J2(Eo + X(k)) (b k % + X - 

k V " 



where bu J>k are creation/anihilation operators of the lo- 
cal mode, and X(k) is the interaction matrix element 
given above (|J). The Hamiltonian Hi oc (||) which de- 
scribes the effective interaction between localized modes 
can be diagonalized using the Bogoliubov transformation 
j3k = u(k)b k + v(k)b 1f The two functions u(k) and v(k) 
are given by: 



a 1 f E +X(k) \ 2 1 f E +X(k) 



- 1 



(6) 



Where E(k), the energy spectrum of the diagonalized 
Hamiltonian is: 



E(k) = ^/E (E Q + 2X(k)) 



(7) 



The ground-state wavefunction of the local-modes is 
given byEj: 



l.„)-n«p(jjV»U) 



(S) 



We show in Fig. 4 X(k) in the (110) direction for a unit 
dipole moment. We now need to fix the size of the dipole 
moment in order to calculate the energy spectrum. 
We would like, according to our definition ot the local 
mode, that the energy cost of flipping the direction on 
a single dipole out of the ground state to be Eq. This 
condition is equivalent to demanding that 2 \X(k = 0)| = 
Eq. We also see from (0) that in order for the dipoles 
to have a gapless mode at k — > we must have that: 
X(k = 0) = —Eq/2. Using this condition, the value of 
Eq = 7K (see previous section) determines the size of 
the dipole-moment as: ~ e • O.OlA. This value is 
indeed smaller than our previous estimation, which was 
an upper bound on the size of dipole moment (|^) . 

From its very definition, the phase modulation in the 
(110) direction of the atomic motion with energy E(k) 
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(Eq.0) should be just the Ti(llO) phonon. In Fig.5 we 
compare the experimental values of Ti(llO) taken from 
neutron scattering with the calculated E(k). The agree- 
ment is excellent for all k. From (Q) and Fig. 4 we find 
that at the edge of the Brillouin zone the energy E(k) 
of the phonon is just the bare energy of the local mode, 
Eq, since X(V2n/a) = and the dipoles have changed 
from the configuration of Fig. 3a to Fig. 3b. We recall that 
the value of Eq that we used was taken from NMR data. 
The agreement between these two independent determi- 
nations of E , that from NMR and that from neutron 
scattering, emphasizes the self-consistency of our descrip- 
tion. We stress that the value of Eq is the only empirical 
input into the calculation, while the functional behav- 
ior is completely defined by the lattice structure and the 
dipolar form of the interactions. 

Our model indicates that only the Ti(llO) phonon 
would be different than obtained from the SCH calcu- 
lation, because it is the only excitation which .-couples 
to the local motion. Indeed, in the experimental this is 
the only phonon branch which is not described well by 
the SCH calculation. The fact that the SCH calculation 
works rather well for other directions, is consistent with 
our picture in which there are no elementary excitations 
of the dipole array in these directions. 



IV. HARMONIC AND LOCAL-MODE 
HYBRIDIZATION 

An equivalent way to describe the mutual influence 
between the local modes and the lattice can be done by 
taking the interaction between the dipoles as resulting 
from an exchange of virtual harmonic transverse fluctua- 
tions of the lattice. An analogous cas&nthat of excitons in 
a dielectric, was treated by HopfieldEj (in that case, the 
interaction is mediated through the exchange of virtual 
photons). The two excitations, the local mode and vir- 
tual harmonic phonon, are then hybridized through the 
same dipolar interaction matrix X{k), which we used in 
the direct interaction picture. In our case, the natural 
choice for the mediating virtual phonon is the Ti(llO) 
phonon, as calculated by the SCH method. The use of 
the phonon calculated by the SCH method is important, 
since this calculation is largely independent of the local 
degrees of freedom described by the local mode. The 
motivation for using this approach is that it allows us to 
obtain an additional branch of the excitation spectrum 
which has observable consequences. 

The Hamiltonian describing the local-mode and SCH- 
phonon hybridization isE3 



H = Hn 



where the three components of (||) are: 
The SCH-phonon 

H Q = y^e(fc)a[a fc 



(9) 



(10) 



where a are Bose operators and e{k) is the SCH- 
phonon spectrum. 
The local-mode 



(ii) 



The part describing the first order (dipolar) cpup un g 
between the phonons and the localized modes is 113 

H c = ( A ( fc ) 6fe + M(fcH) K f + a-k) (12) 
fe 

- (A(fc)fe fc t - fi(k)a k ^ (a k + a f _ fc ) 

where, in the dipolar approximation (and a cubic lat- 
tice) the two functions A and /i are given by A(fc) = 

lE o (-ll$r) * and M*) = -Eo^£$, with X(k) the 
dipole matrix element (Q). This is just the standard 
coupling Hamiltonian of an atom to a transverse pho- 
ton field, which is replaced here by the Ti(110) SCH 
phonon. The total Hamiltonian H does not involve quar- 
tic terms and .©an be diagonalized using the canonical 
transformationtS 



ai 

a 2 = Bat 



Aa k + Bb k + Cat- 



(13) 



where these operators describe the two branches of 
the hybridized energy spectrum. The transformation 
functions— A(k), B{k), C(k), D(k) can be written down 
explicitlycS. The corresponding dispersion relation is 



e 2 (fc) 
E 2 {k) 



= 1 - 



G 



X(k) 



Eq 



1- 



E(k) 
E 



(14) 



which describes two energy branches E\ and Ei- The 
equivalence of this and the treatment in the previous sec- 
tion is due to the use of the same dipole interaction ma- 
trix X(k) in both. With X(k), we can solve (|l|) to find 
the two energy branches: 



E 1 = e{k)/2 
E2 = 2Eq 



(15) 



We see that the energy of the lower branch, E\ is half 
of that of the SCH-phonon for all k. Comparing the SCH 
calculation and the experimental results for the Ti(110) 
phonon, we find that there is indeed a constant, k in- 
dependent ratio between them. This ratio turns out to 
be ~ 1.7 for all momenta, close to the predicted ratio of 
2 ([[5]) (Fig. 6). This small discrepancy is within the un- 
certainty of the SCH calculations and the experimental 
data. Despite this small discrepancy we find the agree- 
ment very satisfying. Our picture of a hybridization of 
a dipole array with the harmonic lattice can therefore 
account for the main effect of the local motion. The up- 
per branch, E 2 , a dispersionless excitation with energy of 
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14K (Fig. 6), was possibly observed in inelastic neutron 
scattering experiments, in which a strong feature was 
observed at an energy transfer of 1.4meV. This feature 
was interpreted by Glydeu in terms of an interference 
effect between phonon modes, which is not inconsistent 
with our hybridization picture. Another possible way to 
observe this excitation would be by Raman scattering. In 
addition, this energy branch is important in the process 
of mass diffusion, and will be discussed in this context in 
a following section. 

We would like to point out that local modes can 
arise also in classical crystals, as a result of a large 
anharmonicitytil. As shown by these authors, these local 
modes can assume many of the roles of vacancies. There 
are several similarities between this work and ours; first, 
the treatment in non-perturbative; second, in a simple 
cubic lattice, they also find that there are two excitation 
branches. These similarities occur despite the different 
basic assumptions of the classical model and the present 
work, which uses the quantum properties of the crystal 
as the starting point. 

V. LOCAL MODES IN HCP 4 HE 

We now turn to discuss the effects arising from the 
local motion in the hep structure. Repeating the calcu- 
lation done in section II, we plot in Fig. 7 the potential 
well along the main directions of the hep crystal, namely 
(1000), (0100) and (0001). We also plot in Fig.7 the low- 
est two energy levels of a one-dimensional Schrodinger's 
equation solved in each of the potential wells. The en- 
ergy difference between these two levels is lowest in the 
(0100) direction, of the order of 16K. Similarly to what 
was stated for the bee structure, the calculation indicates 
the direction along which there will be large amplitude 
local motion. 

As in the bec phase, evidence for the existence of a 
"local mode" comes from NMR measurementsBEj of the 
linewidth (I/T2) with energy ~ 14K. We propose there- 
fore that the NMR experiments measure the energy of 
the bare local mode. 

In contrast to the bec phase, we do not expect the lo- 
cal electric-dipoles in the hep phase to have long-range 
order. This is due to the fact that there is geometric frus- 
tration against an 'antiferroelectric' order in a triangu- 
lar lattice. We calculated the dipolar interaction energy 
for several simple dipole arrangements preserving the net 
zero dipole moment, and did not find any arrangement in 
which this energy was negative. Comparison of the SCH 
calculation of the phonons spectrum with experiment re- 
veals that there is an overall good agreement, with no 
exceptions, such as found for the T^llO) phonon in the 
bec phase. Thus, we conclude that there is no long range 
order within the dipole array like in the bec phase, and 
the local modes in the hep solid are largely uncorrelated. 



VI. ATOMIC SELF-DIFFUSION 

We have seen that in the first excited state of the 
atom in its potential well, the atom oscillates with a 
comparatively large amplitude (~ lA). This motion 
can have an effect on the atomic exchange rate which is 
governed by the overlap of the wave functions of neigh- 
bouring atoms. At low temperature, the atomic ex- 
change, or self diffusion, occurs by tunneling and is small 
( Do ~10~ 9 cm 2 / sec in the hep phasea). Thermally ac- 
tivated diffusion is traditionally attributed to vacancies. 
We would like to show that thermally excited local modes 
can produce self diffusion at rates which compare favor- 
ably with experimental data. 

The theoretical calculation of the zero-temperature ex- 
change tate of atoms in solid Helium is a long-standing 
problemcd. It is extremely difficult to perform accurate 
theoretical calculations of this effect due to the smallness 
of the exchange frequency as compared with the Debye 
frequency. In these calculations, the atoms participating 
in the exchange are treated as tunneling along a one- 
dimensional closed path. At high densities, it is found 
that the effective potential barrier for tunneling is mainly 
due to the potential energy associated with the elastic de- 
formation of the crystal during the atomic exchange. At 
lower densities it is more difficult to calculate the height 
of the tunneling barrier, but empirically it turns out to 
be of the order of the kinetic energy of the atoms (30K 
to 40K). Thus, in order to obtain a rough estimate of the 
exchange rate, it is sufficient to view the problem as that 
of a one-dimensional tunneling of free particles through a 
square barrier. Since we're dealing with low density solid, 
the height of the barrier should be comparable with the 
kinetic energy of the atoms. We would like to compare 
the exchange rate of atoms in the ground state with that 
of atoms in the excited state of the local mode. 

We begin with the hep phase, looking at the (0100) 
direction, that of the largest amplitude atomic motion. 
We take the atom in the ground-state with it's kinetic 
zero-point energy (E^in — 10K) and adjust the bar- 
rier height to reproduce the experimental rate of ex- 
change at low temperatures. We take the width of the 
barrier to be a =3.6A, i.e. the nearest neighbour dis- 
tance. We find that a barrier height of 44K gives a 
transmission probability T ~ 10~ 7 , and a self diffusion 
coefficient ~ (E km /K)a 2 T -10" 9 cm 2 / sec , which 
is consistent with the experimental results. The bar- 
rier height found in this estimation is consistent with 
the empirical value, of the order of the total kinetic en- 
ergy of the atoms in the solid. We now turn to ther- 
mally activated self diffusion. In the first excited state 
of atoms in the well, the total energy of the atom is 24K 
(10K from the ground state + 14K for the first excited 
level of the local mode). The barrier height remains 
unchanged, but the atom can now tunnel from one of 
the lobes of the excited-state wavefunction (see Fig. 2), 
which are approximately lA from its equilibrium posi- 
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tion. Since it tunnels into an identical excited state, the 
effective width of the barrier is now a = 3.6-2=1.61. We 
thus find a transmission probability in the excited state 
r ~ 10~ 2 , and a prefactor of the diffusion coefficient 
D|f ~ (Eki n /K)a 2 T ^10~ 4 cm 2 / sec. This means in our 
model that as atoms are thermally excited out of the 
ground state, the diffusion coefficient would increase ex- 
ponentially with temperature as exp(— Eq / kT) . The 
experimental valuesa are « 10~ 4 cm 2 /sec, and that 
of the activation energy Eq=\'&.§K : which is very close 
to the value of E =14K we took from NMR as describ- 
ing the energy of the local mode. We therefore conclude 
that exchange of atoms occupying a thermally activated 
local mode can account for self diffusion at a rate usually 
attributed to vacancies in the hep phase. 

The rate of self diffusion in the bec phase is an order of 
magnitude larger than in the hep phasefl. In this phase, 
the direction of largest local atomic motion is (100). Re- 
peating the above calculation of the tunneling rate for 
this case, we have the ground-state kinetic zero-point en- 
ergy similar to the hep {Ek in ~ 10K), and barrier width 
along the (100) direction of a =4.12A. For a tunneling 
rate at low temperatures which would be higher by an or- 
der of magnitude compared with the hep phase we need a 
barrier height of ~ 30K, which is again close to the total 
kinetic energy of the atom in the ground-state (~ 34K)El. 
In order to look at thermally activated self diffusion, the 
model used above for the hep phase is inapplicable, since 
the energy of the excited state E 2 = 2E ~ 14°K ( |l5| ) is 
associated with more than one atom. In order to under- 
stand the physical nature of this excitation we plot rib(r) 
(Fig. 8), the spatial extent of the density of atoms excited 
by this branch of the local-modes in the (110) direction. 
This is the Fourier transform of the fc-space density of 
the localized-modes in the upper branch (p"3|): 

n b (k) = (b{b k ) 2 = \C(k)\ 2 (16) 

What is seen in Fig. 8 is that the function rib(r) extends 
over two unit cells in real space, in the (110) direction. 
Since the energy is twice the energy of the bare local- 
mode, it can be interpreted qualitatively as two atoms 
excited in adjacent unit cells, each having an energy 
Eq =7K. This corresponds to a "flip" of the dipole mo- 
ment of these two atoms. In Fig. 9 we show that with this 
excitation, there are 4 atoms oscillating in such a way, as 
to reduce the potential barrier for the exchange of atoms 
1 and 2. This type of atomic exchange resembles phonon 
assisted tunneling, where the correlated " breathing" mo- 
tion of the 4 atoms is locally equivalent to a phonon exci- 
tation at the edge of the Brillouin zone. This type of self 
diffusion was considered in the pastEj, and was recently 
found to be consistent with experiment in bec 4 HeU as 
the dominant majss transport channel, with an activation 
energy of 14.8K&& Thus, altough quantitative calcula- 
tions of the rate of difusion are outside the scope of this 
paper, both the mechanism and the activation energy are 
in accord with experiment. 



VII. CONTRIBUTION OF LOCAL MODES TO 
THE SPECIFIC HEAT 

One of the controversial issues with vacancies in solid 
4 He is that based on the measured values of the activation 
energy, their estimated contribution t^ttie specific heat 
is comparable to that of the phononalSlS. Yet, there in 
no experimental evidence for it. Let us consider the con- 
tribution of the local modes presented here to the specific 
heat. 

Since the energy of the local mode is different in each 
direction of the lattice, the highest possible contribution 
to the specific heat would come from the excitation of 
local modes along the directions where the excitation en- 
ergy is lowest. Because of the restricted solid angle in 
which these modes are active, the fraction of phase space 
occupied by them would be correspondingly small, and 
their contribution to the specific heat would be much re- 
duced compared with that obtained previously. In these 
previous estimations, the point defect was assumed to 
occupy phase space uniformly, thus yielding a large con- 
tribution. To obtain a quantitative estimate of the above 
effect we calculate the ratio between the contribution to 
the specific heat from directional local modes andJiic to- 
tal experimental specific heat of the beco and hepa solid 
phases. We used the energies calculated above for the ex- 
citation energies along the principal directions of the bec 
and hep lattice, and a simple linear interpolation for the 
excitation energies in the intermediate orientations. This 
ratio is plotted against temperature in Fig. 10. For the 
hep phase, the maximum contribution is less than 1%, 
while the bec phase it is ~ 5%. The decrease of this ra- 
tio near T=1.65K is due.-tip the premelting enhancement 
of the total specific heato. These values fall within the 
error bars of the specific heat data, and explain why no 
thermally activated contribution is seen in solid 4 He. 

VIII. CONCLUSION 

In this work we have proposed a new approach to treat 
the local behaviour of 4 He atoms in the bec and hep solid 
phases. We treat the excitations of atoms inside their 
potential well as local modes. The anisotropy of the 
potential renders these modes highly directional. Due 
to the symmetry of the bec phase we propose that the 
local-mode is hybridized with the harmonic density fluc- 
tuations (SCH). The hybridization is described by the 
dipole-dipole interaction and the spectrum of the hy- 
bridized Ti(110) phonon is calculated. An additional 
excitation branch is identified and it is this branch which 
seems to control the anomalously large self-diffusion in 
the bec solid. In the hep phase, the symmetry does not 
allow for correlations of these local modes. Consequently, 
there is no hybridization with the phonons, and the ther- 
mally activated self diffusion in this phase is controlled 
by the energy of the bare local mode. 
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The directionality of the local modes means that their 
contribution to the specific-heat of the solid is negligible. 
We therefore demonstrated that the local mode approach 
can describe experimental data coming from neutron- 
scattering, NMR and diffusion experiments within one 
physical model, while at the same time resolving a long- 
standing discrepancy concerning the specific heat con- 
tribution of point defects. The physical picture of our 
proposed local mode is very different from the classical 
picture of a vacancy. We therefore see no reason to con- 
sider point defects in solid 4 Hc as vacancies, as they can 
be consistently treated as natural excitations of the solid 
lattice, without having to physically remove atoms from 
it. 
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(a) 



-10 12 

o 

Distance from equilibrium (A) 

FIG. 1. The potential- well of an atom in bcc 4 He along dif- 
ferent directions. The energy difference Eo between the low- 
est two energy levels (dashed lines) are: (111)- 59. 5K, (110)- 
27.6K, (100)- 10.6K. 






FIG. 3. The two 'anti-ferroelectric' arrangements of the 
dipoles lying along the major axes of the bcc phase. The 
atoms having the same dipole moments have the same shade. 
The sum of the dipole-dipole interaction (Eq.3) for a unit 
dipole- mo ment are: (a) -0.08, (b) 0.0 (A~ 3 ). 



Distance along (1 00) direction A 
FIG. 2. The probability distribution in the (100) potential 
well (shown in the lower part) . The groundstate wavefunction 
is the solid line while the dashed line is the first excited state. 
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FIG. 5. The experimental data [3] (solid circles) for the 
Ti(llO) phonon compared with the calculation (Eq.7) (solid 
line). Also shown is the energy of the bare local-mode 
E Q = 7K. 
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FIG. 4. The calculated interaction matrix X(k) (Eq.4) as 
a function of the wavevector k, for the three phonon modes 
that could affect a dipolar array. The dipole moment has been 
normalized to give a gapless mode: X(k = 0) = —Eo/2. The 
unit-cell dimensions (A): a = 4.12/2, a' = a\/2. 




0.5 

Wavevector A 

FIG. 6. The SCH calculation [1] (dashed line) compared 
with the experimental results [3] (solid circles) for the Ti(110) 
phonon. Solid line is our calculated spectrum (Eq.7). Also 
shown is the prediction of the hybridization model (Eq.15) for 
the high-energy branch E2 = 2Eq. 
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FIG. 7. The potential-well of an atom in hep 4 He along 
different directions. The energy difference Eo between the 
lowest two energy levels (dashed lines) are: (0100)- 16. OK, 
(1000)- 2I.0K, (0001)- 23.0K. 




Distance along the (110) direction (A) 
FIG. 8. The normalized Fourier transform of the fc-space 
density of the local-mode occupation in the high-energy 
branch E2 (Eq.16), along the (110) direction. 



FIG. 9. Schematic picture of an excitation of the high en- 
ergy branch 2Eo, as two adjacent local modes along the (110) 
direction. The filled circles represent atoms that execute the 
breathing motion, allowing atoms 1 and 2 (open circles) to 
exchange places more easily. 
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FIG. 10. The ratio of the contribution of the local mode to 
the total specific heat in the hep and bec phases [24,25]. For 
clarity, the temperature scale for the bec phase was expanded. 
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